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ABSTRACT 

A maximum likelihood analysis of the NGC 4258 maser positions and velocities reveals a ~ 2a deviation 
from Keplerian motion in the projected rotation curve of the high-velocity features, corresponding to a 
~ 9kms _1 , or 0.8%, flattening of the LOS velocities with respect to Keplerian motion over the range 
of the high-velocity masers. While there are a number of potential explanations for this flattening, we 
argue for pure Keplerian rotation in an inclination-warped disk based on the ability of this model to 
explain a number of otherwise puzzling features of the system. A radial gradient in the disk inclination 
of 0.034 mas -1 is not only consistent with the observed rotation curve, but it generates a bowl along the 
near edge of the disk that naturally explains the otherwise puzzling narrow spread in the declinations of 
the systemic masers. It also explains the existence and location of an apparently recurring flare amongst 
the systemic masers. There is no significant evidence for non-Keplerian rotation in the inclination- 
warped disk. An additional implication of the inclination warp is that the disk rises in front of and 
obscures the central engine at a disk radius of about 8.3 mas, or 0.29 pc. By comparing the observed 
X-ray column to conditions in the disk at this radius, we argue the disk must be atomic at 0.29 pc. 
Hence we conclude that the molecular-to-atomic transition occurs just beyond the outermost maser at 
0.28 pc and from this we infer an accretion rate of ~ 10 aMoyr -1 , where a (< 1) is the standard 
dimcnsionless parameterization of the kinematic viscosity. Our model suggests that most of the observed 
X-ray column arises in the warped accretion disk at 0.29 pc, and that the maser emission is truncated at 
large radii predominantly as a result of the molecular-to-atomic phase transition originally proposed by 
Neufeld & Maloney (1995). The inferred accretion rate is consistent with the jet-dominated accretion 
models of Yuan et al. (2002). 

Subject headings: accretion, accretion disks — galaxies: active — galaxies: nuclei — galaxies: 
individual (NGC 4258) — masers 



1. INTRODUCTION 

NGC 4258 (M106) is a prominent barred spiral galaxy 
(SABbc) located at a distance of 7.2 Mpc (Herrnstein et al. 
1998b). It is one of 12 galaxies included in Seyfert's initial 
catalog of active galaxies (Seyfert 1943) and it has been 
studied extensively over a broad range of frequencies and 
resolutions. On large scales, NGC 4258 is dominated by a 
pair of anomalous spiral arms positioned between the stel- 
lar arms and symmetric with respect to the nucleus. The 
nonthermal emission from these arms spans many decades 
in frequency and is probably energized in shocks formed at 
the interface between the ISM and matter ejected from the 
nucleus (van der Kruit, Oort, & Mathewson 1972; Ford et 
al. 1986). Higher resolution Ha and VLA radio contin- 
uum images reveal that the anomalous arms are fed by 
north-south, multi-strand, twisted jets arising in the ac- 
tive nucleus of NGC 4258 (Cecil, Wilson, & Tully 1992; 
Cecil, Wilson, & De Pree 1995). The jets bifurcate and 
bend sharply about 2.6 kpc from the central engine, pre- 
sumably as a result of interactions with especially dense 
regions in the ISM (Martin et al. 1989). 

Further evidence for nuclear activity comes from IR, op- 
tical, and X-ray observations. Chary & Becklin (1997) de- 
tect a 4.5 mJy compact central source at 2.2 /mi whose 
color is consistent with that of an extincted nonthermal 
source. The near IR emission is unresolved by the 0".6 
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(21 pc) beam. Chary et al. (2000) present higher reso- 
lution images between 1 and 18 /xm, and report nuclear 
emission unresolved by the 0".2 beam with a spectral in- 
dex, s, of 1.4 ± 0.1 (f v — v~ s ) and a 1-20 /jm luminosity 
of ~ 8 x 10 41 ergs -1 . At optical wavelengths, Wilkes et al. 
(1995) detect faint blue continuum emission and broad- 
ened (~ 2500 km s -1 , FWHM) emission in polarized light. 
This is presumably the scattered emission from an ob- 
scured Seyfert 1 nucleus, and NGC 4258 is characterized 
as a hidden broad line region (HBLR) Seyfert 2 AGN. 
Wilkes et al. infer a spectral index of 1.1 ± 0.2 and an op- 
tical luminosity of 5.5 x 10 42 ergs -1 or 5.5 x 10 40 ergs -1 , 
depending on whether the scattering screen is composed of 
electrons or dust. Finally, NGC 4258 has been observed by 
several X-ray telescopes over the last decade (see Table 2 
of Fruscione et al. (2005) for a summary). The most re- 
cent XMM-Newton data of Fruscione et al. (2005) suggest 
a power law spectrum with a neutral hydrogen absorb- 
ing column of 8.6 — 13.2 x 10 22 cm -2 and an absorption- 
corrected 2-10 keV luminosity of 5.4 — 9.1 x 10 40 ergs -1 . 
Taken together, these data imply a bolometric luminosity 
{L bol ) of ~ lO^ergs- 1 . 

While informative, none of these data actually resolve 
the NGC 4258 central engine on sub-parsec-scales, where 
the details of the accretion processes powering the AGN 
are likely to be uncovered. Fortunately, NGC 4258 con- 
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tains a remarkable nuclear water maser, and VLBI ob- 
servations of these dense molecular gas tracers provide a 
unique window into the sub-parsec-scale structure of the 
AGN. The NGC 4258 water maser was initially discovered 
by Claussen, Heiligman, & Lo (1984) and was shown to 
be confined to the inner 3 pc of the galaxy shortly there- 
after (Claussen & Lo 1986). The unexpected discovery 
of maser emission offset ±1000 km s -1 on cither side of 
the brighter emission at the galaxy systemic velocity of 
470 km s -1 (Nakai, Inoue, & Miyoshi 1993) suggested the 
possibility of a rapidly rotating structure. Greenhill et al. 
(1995a) showed that the systemic masers lie along a nar- 
row, east-west structure spanning 4 mas and reported a 
striking linear relationship between position and line-of- 
sight (LOS) velocity that, in the context of disk models, 
suggested either solid body rotation or confinement to a 
narrow range of radii. In parallel with these developments, 
two teams were completing comprehensive, 10-year, single- 
dish monitoring programs of the systemic maser emission. 
Both Haschick, Baan, & Peng (1994) and Greenhill et al. 
(1995b) reported a secular drift in the systemic maser LOS 
velocities of about 9kms _1 yr _1 . In addition, Greenhill et 
al. made the preliminary conclusion that the high-velocity 
masers show no acceleration. Watson & Wallin (1994) 
were the first to propose a self-consistent disk model in 
which the systemic masers lie along a narrow annulus at 
the near edge of an edge-on disk and the high-velocity 
masers are near the diameter defined by the intersection 
of the disk with the plane of the sky (the midline) . 

The VLBI image produced by Miyoshi et al. (1995) pro- 
vides the first unequivocal demonstration that the masers 
trace a well-organized rotating disk (see Figure 1). The 
bulk of high- velocity masers extends from 0.17 to 0.28 
pc on cither side of the systemic masers and apparently 
obeys a Keplerian rotation law. (Note that an isolated 
feature has been detected recently at 0.12 pc (Modjaz et 
al. 2005).) The rotation curve requires a central bind- 
ing mass of 3.8 x 10 7 M Q and a central density of at least 
4 x 10 9 M Q pc~ 3 . This central mass, together with the sys- 
temic maser position- velocity gradient and LOS accelera- 
tions, indicates the systemic masers occupy a narrow annu- 
lus 0.14 pc from the central engine and spanning 8° in disk 
azimuth along the near edge of the disk. At 1100 km s -1 it 
takes individual masers about 10 years to traverse the en- 
tire envelope of systemic emission. The systemic masers lie 
about 0.02 pc south of the disk center as determined from 
the high- velocity masers, requiring that the disk be tipped 
down about 8° from edge on. The lack of structure in the 
vertical positions of the systemic masers implies that the 
aspect ratio of the disk is less than about 0.2% and that 
the portion of the disk traced by the masers is extremely 
thin (Moran et al. 1995). Because the high- velocity masers 
are confined to within several degrees of the midline by 
their small observed accelerations (Bragg et al. 2000), the 
systematic curvature in the declinations of these features 
implies that the disk is warped. The data are well fit 
by a gradual 8° change in the position angle of the disk 
between 0.17 and 0.28 pc (Herrnstein 1997). Herrnstein, 
Greenhill, & Moran (1996) use this warp to explain the 
relative faintness of the blue-shifted emission. An ongo- 
ing multi-epoch VLBA monitoring program has detected 
<~ 30/iasyr -1 proper motions of the systemic masers as 



they sweep in front of the background jet at 1100 kms -1 . 
In addition to confirming the disk geometry, these proper 
motions, together with the accelerations, provide a geo- 
metric distance to the galaxy (Herrnstein et al. 1998b). 

High-frequency VLBA continuum observations have un- 
covered a compact jet elongated along the rotation axis of 
the disk, and well-aligned with the large scale jets at the 
base of the anomalous arms (Herrnstein et al. 1997). The 
~ 0.01 pc offset between the jet emission and the cen- 
tral engine is consistent with optical depth effects along 
a mildly relativistic jet (Blandford & Konigl 1979). The 
southern jet is 5-10 times weaker than the northern jet, 
and this is probably due to thermal absorption in the 
disk. Furthermore, while the high- velocity masers presum- 
ably amplify their own spontaneous emission (the velocity 
gain paths being longest along the midline), the systemic 
masers most likely amplify the southern jet, which lies 
directly behind them. VLBA observations provide a 3-c 
upper limit of 220 /iJy on any 22 GHz continuum emis- 
sion coincident with the central engine (Herrnstein et al. 
1998a). 

Bragg et al. (2000) report the detection with the VLA 
of a flare in the systemic spectrum in 1996, at a velocity 
of 492.6 km s" 1 LSR. The flux density of this flare was 20 
Jy, a factor of 4-5 brighter than the typical peak in the 
systemic emission. The flare maximum was observed on 
29 March 1996. There is some evidence for above-normal 
emission at the same velocity in the preceding observa- 
tion of February 26. However, the observation following 
the flare peak, on 10 May 1996, shows no evidence of the 
flare. From these data, an upper limit on the flare lifetime 
of ~ 75 days is inferred. 

While the warped-disk model derived from the posi- 
tions, LOS velocities, motions, and accelerations of the 
NGC 4258 masers is compelling in its completeness and 
self-consistency, several important questions remain unan- 
swered. In particular, in spite of all that is known about 
this AGN central engine, the accretion rate, M, through 
the disk is still a matter of speculation. The central mass 
together with the estimated bolomctric luminosity indicate 
a highly sub-Eddington system, Lb i ~ 10~ 3 — IQ~ 4 Le, 
the Eddington luminosity corresponding to the central 
mass. In fact, NGC 4258 is quite characteristic of the 
ubiquitous class of relatively faint AGN known collectively 
as low luminosity AGN (LLAGN). A determination of M 
might provide a means to discriminate amongst the various 
accretion modes that have been proposed in NGC 4258, 
and shed some light on accretion processes in LLAGN in 
general. Unfortunately, the infall motions of the gas can- 
not be measured directly in NGC 4258, and estimates for 
M have been indirect and insufficiently precise to constrain 
the accretion models. 

A second category of unanswered questions concerns the 
geometry of the disk itself. The warping in the position 
angle of the disk is immediately obvious because the disk 
mid-line is well-sampled by the high-velocity masers. A 
corresponding warping in the inclination (i) of the disk 
(i.e. i = i(r)) is very likely, as the alternative would im- 
ply that we observe the disk from a very special perspec- 
tive. Unfortunately, the sparse sampling of the disk means 
that observations are not terribly sensitive to an inclina- 
tion warp (Herrnstein, Greenhill, & Moran 1996). There 
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are also questions concerning the precise locations of the 
masers in the disk. It has long been appreciated that the 
high- velocity masers probably occupy the midline because 
this is where the velocity gain paths are longest. Further- 
more, Neufeld & Maloney (1995) explain the outer edge of 
the maser emission in the context of a molecular-to-atomic 
phase transition. However, there is no compelling explana- 
tion for why the systemic masers are so narrowly clustered 
in declination and apparently in radius. Nor is it under- 
stood why the high- velocity features possess an apparently 
abrupt inner edge at 0.17pc on both the red-shifted and 
blue-shifted sides of the disk. 

In this paper we use the shape of the high- velocity maser 
rotation curve to explore the three-dimensional geometry 
of the disk between 0.14 and 0.28 pc. We find that an 
inclination warp consistent with the rotation curve also 
naturally explains the morphology of the systemic masers. 
In addition, we report the discovery of a second flare at 
the same velocity as the first, and show how the inclina- 
tion warp can explain this apparently recurring feature of 
the system. An additional implication of the inclination 
warp is that the disk itself obscures the central engine X- 
ray emission. By comparing the observed X-ray column 
to the predicted disk column, we are able to place tight 
constraints on the accretion rate through the maser disk. 

Section 2 discusses the data and observations utilized 
in the analysis. Section 3 presents the rotation curve fits. 
Section 4 demonstrates that the inclination-warp model 
provides a compelling explanation for the otherwise puz- 
zling geometry of the systemic masers. In Section 5 we 
consider various implications of the inclination warp, fo- 
cusing on the X-ray column and the accretion rate. Sec- 
tion 6 presents concluding remarks. 

2. OBSERVATIONS AND DATA 

The purpose of this analysis is to explore deviations from 
pure Keplerian rotation in the projected rotation curve 
of the high-velocity masers. Isolating such deviations re- 
quires careful statistical analysis, since the precise shape 
of the rotation curve correlates strongly with several ad- 
ditional fitted parameters, most notably the position of 
the disk center. We attempt to minimize statistical uncer- 
tainties in our results by combining a number of datasets, 
listed in Table 1. The first four epochs (designated 'L' in 
the table) were designed to deliver synthesis maps over the 
entire range of maser emission. The principle objective of 
the 'C epochs, on the other hand, was high-sensitivity, 
/zas relative-position images of the central continuum, and 
they yielded synthesis maps for the systemic masers only. 
Each of the epochs included the entire VLBA 4 plus the 
phased VLA acting as a single large aperture. The NRAO 
140-foot and Effelsbcrg 100-meter were included in sev- 
eral of the experiments. The relatively narrow bandwidth 
provided by the VLA complicated the observations consid- 
erably. While the masers span about 150 MHz, the VLA 
is limited to a 50 MHz bandwidth in each polarization. 
Thus, in each observation one polarization was centered 
on the systemic masers while the other polarization al- 
ternated between either the red-shifted and blue-shifted 
masers ('L'), or bands devoid of maser emission, well out- 

4 The National Radio Astronomy Observatory is a facility of the 
Associated Universities, Inc. 



side of the highest-velocity masers ('C'). In the 'L' con- 
figuration we observed a total of four 8-MHz bands, and 
correlated 512 spectral channels across each band giving a 
spectral resolution of about 0.22 km s -1 . The 'C configu- 
ration utilized eight 8-MHz bands and 256 channels, giv- 
ing a spectral resolution of 0.44kms _1 across the systemic 
maser emission. In both frequency setups, the interferom- 
eter was stabilized by applying the phase of a single strong 
systemic maser to the entire spectrum. While this results 
in the loss of all absolute position information, it provides 
relative positions of a small fraction of the 500 ^as beam. 

Our need to phase reference across the full 150 MHz 
spectrum, and across different polarizations, introduced 
a number of complications into the data reduction (Hcr- 
rnstein 1997). Errors, S0 r , in the absolute position of the 
reference maser lead to second order errors in maser po- 
sitions of [Au)/u)]S6 r , where Aw is the offset in frequency 
between the reference and imaged maser, and u> is the ob- 
serving frequency. This amounts to 3.5/xasmas" 1 at 22 
GHz, across the 75 MHz spread of the NGC 4258 maser. 
In this analysis, we reference all the data to the systemic 
maser at 510 km s -1 using the position derived by Her- 
rnstein (1997) via a fringe rate analysis on four VLBA 
epochs: 

a = I2 h 18 m 57.5046 s ± 0.0003 (J2000) 
5 = 47° 18' 14.303" ± 0.003 (J2000) 

The ~ 5 mas uncertainty in the reference maser corre- 
sponds to a negligible 18 [i&s uncertainty across the full 
extent of NGC 4258. In general, broad bandwidth phase 
referencing leads to a heightened sensitivity to all forms 
of delay-like errors. As a result, delay calibrators were 
observed at 10-20 minute intervals throughout the exper- 
iments. These calibrators were also used to calibrate the 
relative phase offsets between the two polarizations, a crit- 
ical step in the calibration process given our need to phase 
reference across polarizations. The VLA was phased-up on 
a nearby, strong calibrator prior to each NGC 4258 scan, 
and typically maintained phase coherence at the 80-90% 
level across the 10-15 minute NGC 4258 scans. The rel- 
ative amount of time on phase-up sources and on NGC 
4258 at the VLA was adjusted in real time according to 
changes in weather conditions and phase stability. 

Each spectral channel was imaged using the standard 
image deconvolution techniques provided in AIPS. Ther- 
mal noises of between 10 and 20 mJy per beam per spectral 
channel were achieved. All features stronger than 2a that 
occurred in two adjacent frequency planes within a beam 
of one another were considered significant. The ccntroided 
positions of these features were determined using a Gaus- 
sian fitting routine in AIPS. Figure 1 includes a spectrum 
derived from the image cubes in this way. The positions 
and velocities of individual masers were derived from this 
spectrum using a crude feature decomposition algorithm. 
Each peak in the interfcromctric spectrum was assumed 
to correspond to a distinct maser feature. The peaks were 
determined from numerical first and second derivatives of 
the spectrum, and the positions and velocities of the result- 
ing maser features are the flux-weighted averages around 
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these peaks. The result is 25-35 systemic masers, 10-20 
red-shifted masers, and 2-4 blue-shifted masers, per epoch. 
Each epoch was referenced to the same point in the sys- 
temic maser emission envelope to ensure consistency. Con- 
tinuum maps were generated by averaging all maser-free 
channels. 

2.1. Single-Dish Monitoring of the Systemic Spectrum 

In hopes of detecting additional systemic flares, the low- 
velocity maser spectrum of NGC4258 was monitored from 
1997 January 27 to 1997 June 24 with the 37-meter diam- 
eter antenna of the Haystack Observatory. Eleven epochs 
were observed. Position-switched total power spectra were 
calibrated using a system temperature referenced to out- 
side the atmosphere and corrected for elevation depen- 
dent antenna gain using a standard model for the an- 
tenna. For each observing track, we observed continuum 
sources 4C39.25 and 1308+326 to calibrate offsets in an- 
tenna pointing. To correct for pointing errors within one 
hour of transit (elevations > 78°), we sampled the peak 
line intensity every ~ 10 minutes and rescaled the spectra 
to achieve a level response. The calibration of the spec- 
trum compiled for each epoch is accurate to <~ 20%. 

We used a 17.8 MHz bandwidth in one circular polar- 
ization and obtained 8192 channel spectra to which we 
applied a boxcar convolving function to obtain a channel 
spacing of 21.7 kHz (0.293kms _1 ). The rms noise after 
removal of a spectral baseline was 0.1 to 0.2 Jy per channel 
for the 11 epochs. 

3. THE HIGH- VELOCITY ROTATION CURVE 

3.1. The Basic Model 



velocities instead of including them in the observed veloc- 
ities. Hence, the relativity-corrected model LOS velocity 
(v' los .) is given by (see, for example, Herrnstein 1997 for 
details) 



Vlost + Vreh 



where 



Vreh — ^ [ V loSi 2 + V roU 2 COS 2 (i) + 2v r oM 



(3) 



(4) 



These model velocities are compared to the observed ve- 
locities (vi) as derived by the classical expression for the 
Doppler shift (optical definition): 

Vi = c— — — , (5) 



where v Q is the rest frequency of the water transition 
(22.23508 GHz), and vi is the observed frequency of the i th 
maser. Our ability to decompose the model velocities into 
classical and relativistic terms (Equation 3) is a convenient 
by-product of using the optical Doppler shift expression. 
The x 2 °f the model fit is given by: 



2 +P 2 Vrot 2 ' 



(6) 



where the summation is over all high-velocity maser fea- 
tures, a Vi are the combined uncertainties in the model and 
observed velocities, and e and (i are geometric correction 
factors discussed below. Because of the excellent spectral 
resolution of the VLBA, a Vi is dominated by uncertainties 
in the model velocities as a result of uncertainties in the 
measured positions of the masers: 



0~Vi — n^roti 

2 Xi - 



(7) 



The shape of the high- velocity rotation curve is explored where a Xi is the uncertainty in the right ascension of the 



using a \ minimization analysis comparing model LOS 
velocities to measured velocities. The details are described 
in this section. In the interest of generality, we revert to 
angular coordinates for the maser positions. Thus, the 
systemic masers are clustered about a disk radius of 3.9 
mas, the high- velocity masers extend from 4.5 to 8.1 mas, 
and at 7.2 Mpc, one mas equals 0.035 pc. The basic thin- 
disk model assumes that the masers obey simple Keplerian 
motion: 

1/2 

(1) 



VroU = 



GM 



yj{xi - x ) 2 + (yi - y a f 



Here, Xi and y, are the angular sky coordinates of the i th 
maser, G is the gravitational constant, M. is the central 
mass (M a ) over the distance (d), and x D and y Q specify the 
disk center. For masers assumed to be along the midline, 
this rotation velocity is related to the LOS velocity (vi OSi ) 
according to 

Vlost = v rot % sin W + Vo, (2) 

where i and v are the inclination (90° corresponding to 
edge-on) and systemic velocity of the disk, respectively. 

Rotation curve models are assessed by comparing vi OSi 
to the observed maser LOS velocities. Special and gen- 
eral relativistic effects play a significant role in the trans- 
formation from observed frequency to LOS velocity for 
the masers. Since the relativistic corrections are model- 
dependent, we choose to remove them from the model LOS 



i th feature, given approximately by O^/^SNR). Here, Q x 
is the x-component of the synthesized beam and SNR is the 
signal-to-noise ratio. We are justified in ignoring the con- 
tribution of a Vi to o~ Vi because of the east- west orientation 
of the disk, e and (3 in Equation 6 correct for the over- 
simplification of assuming that all the masers lie precisely 
along the midline. Observations of high- velocity maser ac- 
celerations indicate a la scatter about the midline of 4.9° 
(Bragg et al. 2000). In the present analysis, we assume 
that the masers are scattered with a uniform probability 
distribution about the midline by an angle ag , taken to be 
independent of radius. In Equation 6, e accounts for the 
fact that this scatter will impose a downward bias on the 
observed LOS velocities. Similarly, j3 specifies the mag- 
nitude of the additional uncertainty in the velocities as a 
result of this scatter. These terms are related to ae as 
follows: 



and 



e~ 1- 



0: 



00 



ere 



1/2 



(8) 



(9) 



The x 2 m Equation 6 is minimized, via a downhill sim- 
plex algorithm, by varying x a , y , v a , M., and i, where 
the disk center is measured with respect to the phase cen- 
ter of the observations, corresponding to the position of 
the maser at 510 kms -1 . Each line in the high- velocity 



5 



spectra was associated with a distinct maser feature, and 
positions were derived as the flux-weighted average across 
each line. The four "L" -setup epochs listed in Table 1 pro- 
vide a total of 99 red features and 17 blue features for the 
fits. Ultimately, it is the relative paucity of blue-shifted 
features that constrains the quality of the fits. 

In practice the high-velocity rotation curve can only be 
used to constrain three of the disk parameters, since M. ' 
and sinz are degenerate, and since the approximate east- 
west orientation of the disk renders the high-velocity fit 
insensitive to y . The results of this section are insensitive 
to y„, and we use y — 0.55 mas as estimated by Her- 
rnstein (1997), based on high-velocity maser spatial sym- 
metry considerations. The inclination can be determined 
from the positions and velocities of the systemic masers. 
Herrnstein (1997) reports i = 82° ± 1 at a radius of 3.9 
mas for a global model incorporating a radius-dependent 
position angle, and a constant inclination. Here, we adopt 
i = 81.4° at 3.9 mas because this results in a consistency 
between the positions of the systemic masers and the bowl 
model derived from the high-velocity data. In fact, the 
conclusions of this section are not substantively affected 
by the choice of inclination over a range of values consis- 
tent with Herrnstein (1997). 

The results of our analyses are shown in Table 2. The 
first two rows of Table 2 show the disk center parameters 
for the best-fitting pure-Keplerian rotation curve model, 
neglecting and incorporating relativistic effects, respec- 
tively. The quoted uncertainties are based on the prob- 
ability density functions of Figure 2, which are derived via 
a Monte Carlo analysis in which x 2 minimizations are per- 
formed on an ensemble of idealized datasets with sampling 
and noise commensurate with the true data. Panels a and 
b of Figure 3 show rotation curves and residuals for these 
models. A ag of 4.8°, corresponding to a velocity scatter of 
~ 4.5 kms -1 at 900 kms -1 , is needed to bring the reduced 
X 2 close to one. Models positioning the masers precisely 
along the midline yield reduced x 2 of over 200. Hence, the 
rotation curve fits provide a highly significant statistical 
detection of the scatter around the midline that is in good 
agreement with the Bragg et al. (2000) measured value of 
4.9°. 

As reported in earlier results, the projected rotation 
curve is at least qualitatively well-fit by a Kcplcrian ro- 
tation law. The data do not yield an unambiguous detec- 
tion of relativistic effects primarily as a result of param- 
eter correlation and the relatively sparse blue coverage. 
These effects are nevertheless incorporated in all subse- 
quent models. 



3.2. Refinements 

We search for deviations from pure Keplerian motion 
by considering four different modifications to the central 
point-mass, Keplerian rotation, constant-inclination disk 
model discussed above. 



3.2.1. Model - Non- Keplerian Rotation 

This first refinement is the most direct and arguably 
the least enlightening probe of the NGC 4258 rotation 
curve. We search for non-Keplcrian rotation across the 



high-velocity masers by replacing Equation 1 by 

Vgm 

VroU = J, (10) 

[(xi - x ) 2 + (yi - y ) 2 f 
and including £ as a fourth parameter in the maximum 
likelihood analysis. The discrepancy between £ and 0.5 
measures the deviation from Kcplcrian rotation in the LOS 
velocities of the masers. Table 2 and panel c of Figure 3 
show the results of such an analysis. The 10% reduction in 
the reduced x 2 is highly significant, and the systematics of 
the blue residuals are noticeably improved. Figure 4 shows 
relative x 2 as a function of £, and Figure 5 shows the £ 
probability density function, derived from the Monte Carlo 
technique described above. We find ~ 2a evidence for de- 
viation from Keplerian rotation ^ 0.5 ) in the NGC 4258 
disk, with a preferred rotation law of y-0.48±0.0i_ This 
amounts to a total flattening in the high-velocity rotation 
curve of about 9.0kms _1 across the high- velocity masers. 

3.2.2. Model 1: - Central Cluster 

Maoz (1995) explored the possibility that the binding 
mass at the center of the NGC 4258 disk is not a su- 
permassive black hole, but a spherically symmetric clus- 
ter of objects. Since such a cluster cannot have abrupt 
edges, Maoz was able to place a lower limit on the cluster 
core density based on the precision of the Keplerian rota- 
tion curve. However, this analysis relied on the relatively 
crude upper limits on the deviation from Keplerian rota- 
tion available at the time, and we duplicate it here using 
the x 2 formalism developed in the previous section. 

Following Maoz, we consider the possibility that all the 
binding mass resides in a spherically symmetric Plummer- 
model cluster at the center of the disk. The mass distri- 
bution is then parameterized in terms of the core density 
and radius (p c and r c ) as follows: 

-3/2 

(11) 



M{< r) = r 



1 + 

Tc 2 



The cluster model is implemented in the x 2 analysis by 
associating M(< r{) with M in Equation 1 and adding 
p c as a fourth parameter in the fit. Once again, Table 2 
indicates a significant improvement in the reduced x 2 ° ver 
pure Keplerian models, and Figure 3 shows the fit and its 
residuals. 

A disk obeying perfect Keplerian rotation would not 
possess a X 2 (Pc) minimum. Instead, the relative % 2 would 
asymptotically approach zero for large values of p c and 
such an analysis would provide only a lower limit on p c . 
The minimum in relative x 2 {Pc) evident in Figure 4 there- 
fore implies a significantly flattened rotation curve. Specif- 
ically, we find p c = 5.7 ± 1.4 x 10 11 M Q pc~ 3 (Figure 5), a 
core radius of 0.78 mas (= 8.4 x 10 16 cm), a total clus- 
ter mass between the inner and outermost masers of 
1.1 x 1O 6 M , and a resulting rotation curve flattening of 
10.9 kms -1 across the high- velocity masers. 

3.2.3. Model 2 - Massive Disk 

We now consider the possibility that mass in the outer 
accretion disk itself is responsible for flattening the rota- 
tion curve. The maser observations provide direct evidence 
that the outer disk is thin, cool, and dense. For a thin, 
isothermal, steady-state accretion disk in hydrostatic equi- 
librium (c.f. Frank, King, & Raine 2002), the scale height 
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is c s r 3 / 2 {GM) 1 / 2 , and the midplane density as a function 
of radial distance r can be written (Neufeld & Maloney 
1995): 

GM M 

Pmld - 3^)V2 ar 3 Cs 3 [LZ > 

where M a is the central mass, M is the accretion rate, c s 
is the sound speed, and a is the Shakura-Sunyaev param- 
eter of the kinematic viscosity. Integration of Equation 12 
gives the mass of disk, Md, within the inner and outer 
boundaries of the maser emission, r\ and r 2 , 



Md = 8.3 x 10 



4 M 



My 2 (rV 2 -r 



V 2 ) 



a 



0i 



(13) 



where r\ and r 2 are in pc, c s is in kms -1 , M is in solar 
masses per year and M Q is in solar masses. The disk mass 
internal to the i th maser feature (M Di ) is estimated from 
Equation 13, using r 2 = Vi and r\ = r min (the radius of 
the innermost maser), and assuming c s = 2.3 kms -1 , cor- 
responding to a temperature of 800 K. Theoretical mod- 
els predict that temperatures between 500 and 1000 K are 
most conducive to the formation of water masers. 

A massive disk is included in the rotation curve model 
by incorporating Ma -1 as a fourth parameter in the fit 
and replacing M. in Equation 1 by M. + M^Jd 1 where 
the distance, d, is taken to be 7.2 Mpc. The reduced x 2 
improvement is comparable to the cluster model results 
(Table 2), and the rotation curve and residuals arc in- 
cluded in Figure 3. As with the central cluster analysis of 
the previous section, the pronounced minimum in relative 
X 2 (Ma _1 ) evident in Figure 4 implies a significantly flat- 
tened rotation curve. The x 2 analysis favors an accretion 
rate of 6.5 ± 1.3 x 10~ 2 aM Q yr" 1 (Figure 5), which leads 
to a total disk mass of 8.9 x 1O 5 M and a ~ 10.1 km s -1 
flattening in the rotation curve across the high-velocity 
mascrs. Note that for this disk mass the midplane densi- 
ties at the inner and outer radii of 0.14 and 0.28 pc are 
4 x 10 11 and 0.8 x 10 11 cm" 3 respectively, about an or- 
der of magnitude above the normally preferred range for 
maser emission of 10 8 to 10 10 cm~ 3 . Hence, the disk con- 
ditions implied by the existence of maser emission make 
this model unlikely. 

3.2.4. Model 3 - Inclination Warp 

The disk geometry is parameterized in terms of two Eu- 
ler angles: the position angle (<j>) is azimuthal about the 
LOS while the inclination angle is polar. A warped disk 
can be described by allowing these angles to vary with ra- 
dius. While the specific functional relationship between i 
and 4> and r will depend on the perspective of the observer, 
it can be expected that both Euler angles will have some 
radial dependence. As discussed in Section 1, the declina- 
tions of the high-velocity masers are well-fit by a position 
angle obeying a second order polynomial in radius. Thus, 
it seems likely that the inclination of the disk will also vary 
with radius. Such an inclination warp will manifest itself 
as an apparent deviation from pure Keplcrian rotation in 
the LOS velocities of the high-velocity mascrs. 

We quantify this effect by postulating the following (ad 
hoc) linear inclination warp model: 

di 



+ ^(ri-ro), 



(14) 



where ii is the inclination of the disk at the radius of the 
i th maser. As discussed above, i D and r a are fixed at 81.4° 
and 3.9 mas, respectively, based on observations of the 
systemic masers. The inclination warp is added to the pro- 
jected rotation curve fit by replacing i by ii in Equation 2 
and including di/dr as a fourth free parameter. This leads 
to a significant improvement in the reduced x 2 over pure 
Keplcrian, constant-inclination models (Table 2) and the 
model and residuals are shown in Figure 3. The x 2 anal- 
ysis prefers an inclination gradient of 0.04 ± 0.01 mas -1 
across the high- velocity masers (Figures 4 and 5). Con- 
sistent with previous results, this corresponds to a rota- 
tion curve flattening of about 9.6 km s -1 across the high- 
velocity emission. In this model, the inclination of the disk 
changes from 81.4° at the radii of the systemic masers 
(3.9 mas), to about 91° at the radius of the outermost 
high- velocity feature (8.1 mas). The magnitude of the in- 
clination warp is approximately equivalent to that of the 
position angle warp. 

We emphasize that, unlike the other models under con- 
sideration, the inclination warp model does not indicate 
non-Keplerian motions in the NGC 4258 disk, but rather 
radius dependent projection effects in the disk. 



3.3. Discussion 

While there is compelling evidence for a <~ 8kms _1 flat- 
tening in the projected rotation curve of the high-velocity 
masers, its physical origin is ambiguous. In the previous 
sections, we have examined three possible mechanisms, 
each in principle capable of generating the observed ro- 
tation curve. Not surprisingly, the parameters defining 
these three mechanisms are highly correlated with one an- 
other. Furthermore, because an inclination warp is capa- 
ble of both flattening and steepening the rotation curve, 
the actual disk could in fact incorporate any combination 
of these effects. This is demonstrated graphically in Fig- 
ure 6, which shows contours of constant x 2 throughout the 
di/dr-Ma^ 1 plane. The correlation between the parame- 
ters is clear, and arbitrarily large Ma -1 can be supported 
by positing negative values for di/dr without incurring too 
dramatic an increase in x 2 - This limits our ability to place 
firm limits on either Ma" 1 or p c from the high- velocity 
rotation curve alone. We note in passing that the overall 
minimum in the joint x 2 of Figure 6 suggests a detection 
of the effects of finite mass in an inclination-warped disk. 
However, the detection is very weak (a Ax 2 of ~ 2 for the 
inclusion of M as a new parameter) and the densities im- 
plied by the preferred accretion rate are inconsistent with 
the existence of masers. 

The three mechanisms considered above are not in- 
tended to be an exhaustive list of the mechanisms capable 
of flattening the projected rotation curve. For example, it 
might be possible to fit the rotation curve without violat- 
ing the acceleration data by systematically displacing the 
innermost masers farther from the midline than the outer 
masers. We merely note in passing that such a geometry 
is counterintuitive in that the envelope of maser emission 
is expected to grow with radius, along with the coherent 
velocity gain paths. 
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4. THE SYSTEMIC MASERS AND THE INCLINATION WARP 

The previous section suggests that additional informa- 
tion is required to discern the true origin for the flattened 
rotation curve. In this section we show that the inclina- 
tion warp derived from the high-velocity masers provides 
a compelling explanation for the morphology of the sys- 
temic masers. Based on this we conclude that inclination 
warping is indeed the source of the apparently flattened ro- 
tation curve. We begin by first demonstrating that there 
are substantial unanswered questions surrounding the ge- 
ometry of the systemic features. We then show that the 
inclination warp answers these questions convincingly. 

4.1. The Systemic Masers 

The initial VLBA observations of the NGC 4258 masers 
(1994 April) failed to detect any vertical structure in the 
systemic masers at the 10/zas level, leading Moran et al. 
(1995) to conclude that the disk aspect ratio (H/r) must 
be less than about 0.25%. We have repeated this analy- 
sis using epochs 6-8 in Table 1. These observations in- 
clude the Effelsberg 100-m antenna and are about three 
times more sensitive to vertical structure. The results are 
shown in Figure 7. The slope of the best-fitting line sug- 
gests the systemic masers trace a disk with a 5.6° position 
angle. The linear fit has a reduced x 2 01 82 when the 
y-coordinate instrumental uncertainties are used, indicat- 
ing a highly significant detection of vertical structure in 
the systemic masers. Specifically, a reduced \ 2 °f one is 
achieved by adding in quadrature an additional scatter in 
y (a' y ) of 8/ias or 0.00024 pc, indicating that H/r < 0.2% 
at the radius of the systemic masers. 

For a constant inclination disk, the narrow spread in y 
of the systemic masers is puzzling in light of the fact that 
the high- velocity masers extend from 4.7 to 8.1 mas. It is 
even more surprising when we consider that the systemic 
masers apparently amplify background continuum emis- 
sion from the compact nuclear jet. Figure 8 shows, for a 
constant-inclination disk, z (distance along the LOS) ver- 
sus y along the radial perpendicular to the disk midline, 
with disk azimuth (9) of 270°. 9 is measured in the direc- 
tion of rotation from the ascending node of the disk, such 
that 6 — for the red-shifted masers. The LOS velocities 
of the systemic masers range from about v Q — 70kms _1 to 
v Q + 70 km s _1 indicating that they occupy a ~ 7° arc cen- 
tered about 9 = 270°. This is a natural place for masers 
to occur since, together with the midline, it is the location 
in the disk where dvi os /dz is a minimum. Furthermore, 
as long as the dependence of position angle on radius is 
not too steep at small radii, jet emission aligned with the 
rotation axis of the disk lies behind the 9 — 270° radial 
irrespective of position angle, and masers located here are 
well-situated to amplify background jet emission. 

The vertical dash in Figure 8 marks the position and y 
spread of the systemic masers. Continuum emission has 
been detected in seven of eight VLBA epochs. The north- 
ern jet emission is variable in both strength, ranging from 
less than 2mJy to 6mJy, and in position, with northern 
offsets between 0.35 and 1.0 mas. The southern jet has 
remained relatively stable at 0.5 mJy and 1.0 mas South 
of the disk center (Hcrrnstcin et al. 1997, 1998a). The 
observed range of northern jet emission is marked by the 
upper shaded oval in Figure 8. The southern shaded oval 



is the reflection of the northern emission through the disk 
center, and the asterisk marks the location of the observed 
southern emission. Herrnstein, Grccnhill, & Moran (1996) 
demonstrate that any portion of the disk directly exposed 
to the central X-rays will develop a thin layer of ionized 
gas with 2 2- GHz free-free optical depths of > 2 that scales 
as <~ r~ 2 . They further propose that the southern jet ap- 
pears weaker and is displaced southward as a result of ther- 
mal absorption in the disk. Figure 8 shows that if the jet 
is symmetric, the brightest 22 GHz emission lies directly 
behind the systemic masers suggesting that these masers 
amplify southern jet emission that is otherwise obscured 
from view by the disk itself. In this scenario, the extremely 
bright systemic masers are only mildly saturated, and their 
variability is naturally explained by variability in the back- 
ground continuum. Both the northern jet continuum and 
the systemic masers vary on timescales of weeks. 

Clearly there is a serious inconsistency between the 
background amplification models and the observed nar- 
row y extent of the systemic masers at 9 ~ 270°. While 
the background 22 GHz continuum radiation extends 
to —1.0 mas with respect to the disk center, the sys- 
temic maser emission is narrowly confined around y = 
—0.55 mas. Furthermore, symmetry with the northern 
jet emission would imply that amplifying systemic masers 
should range in y from -0.35 mas to -1.0, about a factor of 
ten larger than the observed spread. Of course, the bright 
systemic masers require both background seed photons and 
an inverted population of water molecules. However, it is 
very difficult to construct a pumping mechanism that leads 
to such a narrow annulus of inverted molecules. Ncufcld 
& Maloney (1995) show that the high-velocity masers are 
pumped by direct exposure to the central X-ray emission, 
as a result of the position angle warp. No plausible incli- 
nation warp would yield such a narrow cusp of pumped 
water on the near edge of the disk. One possibility is that 
the outer cool disk is truncated within 3.9 mas, and that 
the systemic masers trace the narrow annulus of pumped 
gas along the inner edge of the disk. However, we know of 
no physical explanation for such a truncation. 

The narrow spread in the systemic declinations can 
also place more quantitative restrictions on a constant- 
inclination disk. For a disk with no inclination warp, the 
observed scatter in y places a firm upper limit on the ra- 
dial scatter (oy) of the systemic masers. For a disk tipped 
down 8° from edge-on, a' y ~ 0.14oy, and the detected ver- 
tical scatter in the systemic masers requires oy < 0.06 mas 
(= 0.002 pc). Thus, for a constant-inclination disk, the 
narrow spread in y of the systemic masers requires a frac- 
tional radial spread of less than 1.5% at a radius of 3.9 mas. 
This requirement can be tested directly using the velocities 
of the systemic masers. 

If the systemic masers are indeed distributed az- 
imuthally in a Keplcrian disk, then 
dvi os GM 1/2 sin(z) cos(0) 

— - — ~ = constant. (15) 

ax r s d/ 
where <j) is the disk position angle. The instrumental un- 
certainties in the velocity are negligible, and we express 
the uncertainty in v, oy, as the quadrature sum of two 
terms: 



oy = 270a a x 2 



3 oy , 
- — (x - 
2 r v 



(16) 
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where a x is the instrumental uncertainty in right ascen- 
sion and x' a — x a + r a sm(<f>) cos(zo) is the projected near 
edge of the disk where the observed LOS velocity is pre- 
cisely the systemic velocity and is independent of ra- 
dius. 270 km s -1 mas" 1 is the well-constrained prior on 
the position- velocity slope based on numerous, past obser- 
vations. The first term in the sum is the straightforward 
propagation of position uncertainties through to a v , while 
the second propagates u r through to a v . 

The fit in the upper panels of Figure 9 uses a r /r — 0.015 
in conformation with the upper limit on the fractional scat- 
ter in radius from the observed scatter in y. The resulting 
reduced \ 2 of 40 demonstrates a significant inconsistency 
in the distributions of the declinations and velocities of the 
systemic masers in the context of a constant-inclination 
disk. In the lower panels a reduced x 2 of one is achieved 
for a r /r = 0.13, a factor of 10 larger than indicated by 
the declination data. Hence, a constant-inclination disk 
model is unable to accommodate both the extremely nar- 
row scatter in declination and the more substantial velocity 
deviations of the systemic masers. 

4.2. The Inclination Warp: Masers in a Bowl 

In the rotation curve analysis of Section 3.2.4 we found 
di/dr — 0.04 ± 0.01 mas -1 . In the sections that follow, we 
will assess this estimate based on a variety of considera- 
tions involving the systemic masers. The resulting refined 
estimate for the inclination warp, di/dr', of 0.034mas _1 
is the best-guess estimate for the warp based on the high- 
velocity and systemic data together. In motivating the 
inclination warp in the present section, we are justified in 
using this value for the warp, as it is well within the un- 
certainties indicated by the high velocity data alone (Fig- 
ure 5). 

Figure 10 shows the 6 = 270° cross-section of a disk 
with di/dr — 0.034mas~ 1 . At the observed y-offset of the 
systemic masers, marked by the vertical dash in the figure, 
9y/9r|(g =2 70) ~ 0- I n other words, for a di/dr consistent 
with the high-velocity rotation curve, the systemic masers 
lie at the bottom of a bowl along the near edge of the disk, 
as viewed from our perspective. That the systemic masers 
lie on the disk at all at y = —0.55 is critically sensitive to 
our choices for y a , i a , and r a , and indeed determined the 
values used for these parameters in the fits of Section 3.2.4. 
However, the fact that these high- velocity fits lead to a ge- 
ometry in which the systemic masers lie at a tangent point 
to the disk is either a remarkable coincidence or of physical 
significance. 

In fact, there is a compelling reason the brightest sys- 
temic emission should occur at the bottom of such a bowl. 
Figure 1 1 is a schematic representation of the bowl cross- 
section. The upper surface of the bowl is directly exposed 
to X-rays from the central engine and will therefore possess 
the aforementioned thin surface layer of ionized gas with 
optical depths > 2. Any molecular gas suitable for sup- 
porting masers will lie underneath this ionized gas. Three 
radii are labeled in the Figure: B marks the bottom of the 
bowl, while points A and C delineate the range of radii 
with unattenuatcd sightlines to both the background con- 
tinuum and us. Water molecules internal to A will have 
an unobstructed view of the background jet. However, any 
mascr emission generated in this gas must traverse two 



ionized surfaces of the bowl along the LOS to us and will 
therefore be attenuated by a factor > 50. Similarly, while 
molecular gas external to C will have an unobstructed LOS 
to us, it will have an equally attenuated view of the back- 
ground continuum and will effectively be limited to self 
amplification. By analogy with the self-amplified high- 
velocity masers, such emission is likely to be at least 20 
times weaker than the background amplified emission in 
the shaded region at the bottom of the bowl. 

A significant shortcoming of the constant-inclination 
disk model is its inability to explain the narrow y- 
confinement of the systemic masers (Figure 8). In con- 
trast, the bowl model predicts that the brightest systemic 
emission ought to lie at the bowl bottom in a narrow dec- 
lination band on the sky, corresponding to the thickness of 
the masing region along the bottom of the bowl (shaded re- 
gion in Figure 11). Figure 10 shows that this is indeed the 
case for the di/dr estimated from the high- velocity data. 
In addition to predicting that the brightest systemic emis- 
sion arises in the bowl bottom, the bowl model also allows 
for much (~ 20 — 50 times) fainter systemic emission at 
larger or smaller radii, as described above. Such emission 
is not ruled out in the present analysis. 

The top panel of Figure 12 once again shows systemic 
maser x versus y. The solid line demarcates the bowl bot- 
tom, and demonstrates in detail the coincidence in decli- 
nation of the systemic masers and the bottom of the bowl. 
The parallel dotted line is intended to represent the ap- 
proximate maser-layer thickness needed to accommodate 
all the systemic emission. It corresponds to a la thickness 
of 15/ias. It is larger than the la thickness inferred from 
a straightforward linear fit to the data (Figure 7) because 
the bowl bottom is not parallel to this best-fitting line. 
Figure 12 suggests that the declinations of the systemic 
masers require a masing layer of ~ 15/ias (la) at the bot- 
tom of the bowl. The bottom panel of Figure 12 shows 
radius as a function of x for the three points identified in 
Figure 11. There is a ~ 10% gradient in the radius of 
the bowl bottom across the envelope of systemic emission. 
The inner and outer edges of the bowl (A and C) are lo- 
cated 0.77 and 1.23 times the bowl bottom radius for a 
maser layer thickness of 15/ias. 

A second serious deficiency of the constant-inclination 
disk is an inconsistency between the extremely narrow dec- 
lination scatter and relatively larger velocity spread of the 
systemic masers (Figure 9). Figure 12 shows that masers 
along the bowl bottom, and therefore confined to a nar- 
row band of declination on the sky, can, in principle, reside 
over a large range of radii in the disk. Because the bowl 
model removes the proportionality between y and r char- 
acteristic of constant-inclination models, it also removes 
any inconsistency between the declinations and velocities 
of the systemic masers. 

This conclusion is quantified in Figure 13, which shows 
actual systemic maser position-velocity data superposed 
on the predicted LOS velocity for the three annuli iden- 
tified in Figure 11, assuming a maser layer thickness of 
15/ias. The precise relationship between the model x 
and v depends on both the inclination and position an- 
gle warps. As in previous works (Herrnstein, Greenhill, & 
Moran 1996), we estimate the latter warp by assuming the 
high- velocity masers lie along the disk midline, and fitting 



9 



a second order polynomial through their positions. Using 
the first four epochs of Table 1, we find the data are well 
fit by 4> = 0.425 - 0.088r + 0.0038r 2 (r in mas), in good 
agreement with past estimates for the position angle warp. 

If the systemic masers were positioned exactly along the 
bowl bottom, we would expect the data to lie along the 
line B in Figure 13. However, given the scatter in the 
systemic declinations and the bowl curvature determined 
from the high-velocity data, the systemic masers could in 
principle lie anywhere in the range of radii shown in Fig- 
ure 12. Hence, the lines A and C show as a function of 
x the envelope of LOS velocity that is consistent with the 
bowl model and the observed scatter in y. Whereas the 
constant inclination disk is unable to handle the structure 
in the systemic x—v data, and in particular the structure 
between —0.05 and —0.25 mas (Figure 9), the bowl model 
has no trouble accommodating the observations. 

The model velocities in Figure 13 were generated using 
v Q = 473.5 km s" 1 and x — —0.19 mas. The convergence 
point of the three velocity curves in the figure marks the 
point along the near edge of the disk where the LOS ve- 
locity gradient across the bowl bottom and along the line 
of sight reaches a minimum of < 2kms -1 . A notewor- 
thy feature of the model is that the LOS velocity gradient 
across the bowl is a minimum for a velocity and a po- 
sition that are ~ 20 kms -1 and <~ 0.1 mas to the red of 
the systemic velocity and projected near edge of the disk 
[x' ~ —0.15 mas), respectively This is also evident in 
Figure 14, which shows lines of constant LOS velocity in 
the bowl assuming x a — —0.19 mas. The iso- velocity lines 
are most nearly vertical for velocities ~ 20 kms -1 red of 
v Q . In the next section we make an association between 
this velocity convergence point and the observed systemic 
flare. 

4.3. The Inclination Warp: The Systemic Flare 

The 1996 systemic maser flare arose at a velocity of 
about 494 kms -1 (Bragg et al. 2000). Figure 15 shows 
the single-dish Haystack detection of a second flare in the 
NGC 4258 systemic spectrum. The flare was first unam- 
biguously evident on day 77 of 1997, at which point it was 
near its maximum flux density of about 17 Jy. It returned 
to a more standard strength of ~ 5 Jy by day 175 of 1997. 
The implied lifetime of this flare is thus < 100 days, some- 
what longer than the 75-day upper limit indicated by the 
Bragg et al. (2000) observations of the 1996 flare. The ccn- 
troided velocity of the flare reddens from 494.1 to 496.4 
kms -1 over the course of the observations, and is well- 
fit by a secular drift of 8.21 ± 0.82 km s -1 yr^ 1 , which is 
consistent with previous detections of accelerations in the 
systemic masers (see Figure 16). Coincidentally, VLB A 
epochs BH25a, b, and c (Table 1) occurred near the maxi- 
mum of the 1997 flare and detected a peak flux density of 
17 Jy and flare velocities between 493.3 and 495.0 kms -1 
(Figure 17). 

Figure 18 shows spectra for both the 1996 and 1997 
flares. The flare flux density maximums occur between 
492.6 km s" 1 (in 1996) and 494.6 kms -1 (at the end of 
the 1997 flare). The line profile of the 1996 flare was in 
fact not well-fit by a single feature. In a three-feature fit, 
the two features significantly brighter than standard sys- 
temic emission are centered about 488.4 km s -1 (14 Jy) 



and 492.6 km s -1 (20 Jy). The impressive agreement in 
the velocities of the two flares suggests that the NGC 4258 
flare is actually a recurring event at a particular point in 
the disk. That is, flaring is observed when a suitable patch 
of enhanced density gas moves through a zone of optimum 
amplification. For x Q — — 0.19mas and v = 473.5kms~ 1 , 
Figure 13 shows that the flare position corresponds with 
the velocity convergence point of the disk 20 kms -1 and 
0.1 mas to the red of disk center. 

Hence, we propose that the NGC 4258 flare is recurring, 
and arises as a result of the bowl-like geometry along the 
front of the disk. At vi os rts v Q + 20 kms -1 the inclina- 
tion and position angle warps conspire to bring the LOS 
velocity gradient along the LOS to < 2kms~ 1 across the 
entire ~ 2 mas width of the bowl. We speculate that flares 
arise as a result of the serendipitous alignment and coher- 
ent self- amplification of masers along the bowl bottom at 
<~ v Q + 20km s -1 . The time scale between flares of about 
one year suggest that the characteristic length scale be- 
tween gas density peaks is about 2 x 10 16 cm. We note 
that Boboltz et al. (1998) have shown that a 1990 flare in 
the W49N water maser is well-modeled by the same type 
of chance alignment proposed here. 

The inclination and position angle warps indicated by 
the positions and velocities of the high- velocity masers un- 
ambiguously predict a velocity convergence point along the 
near edge of the disk. However, the preferred x and v a of 
the high- velocity fits do not lead to a precise coincidence 
between the velocity convergence point and the location 
of the flare. The final entry of Table 2, labeled di/dr', 
is the best-fitting inclination-warp model that also leads 
to center coordinates consistent with the flare. It is ar- 
rived at using the following iterative scheme. The high- 
velocity fitting is first performed with no constraints on 
x a or v , yielding the solution di/dr. The velocity con- 
vergence point, which depends on the precise shape of the 
warp, the central mass, and the disk center coordinates, is 
computed for this model. x{ and v£, the disk center co- 
ordinates that locate the flare at the velocity convergence 
point while holding all other model parameters fixed, are 
estimated, and the high-velocity fit is rerun around x£ 
and v[ priors. The sequence is repeated with increasingly 
narrow prior distributions about x{ and v{ until conver- 
gence is achieved. In this case, convergence indicates the 
existence of an x^ and vl for which the best-fitting cen- 
tral mass and warp model (as determined from the high- 
velocity masers) also leads to an alignment between the 
velocity convergence point and the flare. 

Table 2 shows that there is a ~ 2cr discrepancy be- 
tween the unconstrained high-velocity fit and the flare- 
constrained fit. There is also a significant degradation in 
the reduced x 2 of the flare-constrained fit. Past efforts at 
global disk fitting have uncovered comparable inconsisten- 
cies between the systemic and high velocity masers, and 
this is summarized in Figure 19. The diamonds denote 
the maximum-likelihood centers for fits considering only 
the high-velocity masers. The squares mark the center 
values for fits incorporating both the systemic and high- 
velocity masers. The open and filled symbols represent 
position-angle warp and inclination- and position-angle 
warp models, respectively. The open square is derived 
from true global fit incorporating all masers, while the 
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filled square corresponds to the di/dr 1 solution discussed 
above, and utilizes the high-velocity masers and the lo- 
cation of the flare. The discrepancy between the open 
diamond and open squares illustrates the tension between 
the high-velocity and systemic masers in the context of 
a position-angle warp. The solid symbols indicate that 
by including an inclination warp and locating the flare at 
the velocity convergence point of the disk, we have signif- 
icantly reduced, but not removed this tension. 

4.4. Discussion 

In this section we have attempted to demonstrate that 
the inclination warp derived from the projected high- 
velocity rotation curve naturally explains a number of oth- 
erwise puzzling features of the systemic masers. In par- 
ticular, we have shown that the bowl geometry provides 
an explanation for the extremely tight declination spread, 
and broader velocity spread of the systemic masers. Fur- 
thermore, we have argued that the combination of cen- 
tral X-ray heating, background amplification, and ther- 
mal absorption makes the bowl bottom an especially at- 
tractive location for the systemic masers to arise. Finally, 
we have speculated that the apparently recurring flare in 
NGC 4258 arises at the point in the bowl where the ve- 
locity gradient along the LOS approaches zero. Based 
on these results, we conclude that di/dr = di/dr' = 
0.034 mas^ 1 and that the inclination warp is in fact pre- 
dominantly responsible for the flattened high-velocity pro- 
jected rotation curve. Figure 20 shows three perspectives 
of this disk. The disk midline and the 9 = 0° and 9 = 270° 
radials are overlaid. The peanut-shaped curve traces the 
location of the disk tangent, as seen from our perspective. 
As prescribed by the model, the systemic masers have been 
distributed along the disk tangent at the near edge of the 
disk, while the high velocity masers have been placed along 
the midline. 

5. IMPLICATIONS 

5.1. The Rotation Curve Revisited 

Of the three specific refinements to a constant- 
inclination, centrally bound disk considered in Section 3.2, 
it was the inclination warp alone that was capable of steep- 
ening the rotation curve and that prevented us from plac- 
ing interesting limits on £, p c , and Ma^ 1 . The dotted 
curves in Figure 4 show x 2 (£); X 2 (Pc), and ^{Ma -1 ) for 
a disk with di/dr = 0.034mas~ 1 . Figure 5 also shows 
the probability distribution for £ in the presence of an 
inclination warp. The maximum-likelihood rotation law 
for a inclination-warped disk is v rot oc r -°- 49 ±o.oi^ an( j. 
there is no significant evidence for deviation from pure- 
Keplerian rotation. The asymptotic behavior of the rela- 
tive p c and M« _1 % 2 further reinforces the conclusion that 
the inclination warp alone fits the projected high-velocity 
rotation curve. For di/dr = 0.034 mas -1 , we find that 
Pc > 10 12 M Q pc- 3 and M < 5 x lO- 2 aM yr" 1 . 

Maoz (1998) has developed a formalism for estimating 
cluster lifetime upper limits {T max ) by considering both 
collision and evaporation timescales for a variety of clus- 
ter constituents. For NGC 4258, he finds T max < 1 Gyr for 
p c > 8 x 10 11 M Q pc~ 3 . Hence, the lower limit on p c de- 
rived from the high-velocity rotation curve indicates that 



if the disk is bound by a central cluster, the cluster life- 
time would have to be considerably shorter than the age 
of the galaxy itself. This strongly favors the presence of a 
supermassive black hole in NGC 4258. The upper limit on 
Ma" 1 is less interesting since it is three times larger than 
what would be expected for advection-dominated accre- 
tion flow models (ADAF; Gammie, Narayan, & Blanford 
1999), which are generally the least efficient, highest M 
solutions. 

5.2. The Disk Geometry Revisited 

For a disk with no position angle warp, the narrow 
spread in the declinations of the systemic masers implies 
that these features are tightly clustered about a radius of 
0.14 pc (3.9 mas). A number of theories have been pro- 
posed to explain the inner cutoff in the systemic emis- 
sion. Neufeld & Maloney (1995) suggest that the inner 
edge of systemic emission might be due to a flattening of 
the disk and a resulting decrease in the X-ray heating in- 
side of this radius. More recently Yuan et al. (2002) have 
fit the NGC 4258 SED with a jet-dominated AGN model 
in which the outer thin disk transitions to an ADAF at 
a radius r tr . The NGC 4258 SED alone does not place 
tight constraints on rt r , and it is tempting to associate rtr 
with the inner edge of systemic emission. However, if the 
inclination warp model is correct then the systemic masers 
are clustered along the bottom of a bowl along the near 
edge of the disk, and there is no additional explanation 
needed for their narrow spread in y. In this model, there 
is nothing distinctive about 0.14pc except that from our 
perspective it corresponds to a tangent point to the disk. 

Similarly, efforts have been made to explain the observed 
range of high-velocity emission from 0.16 to 0.28 pc (4.5- 
8.1 mas) on either side of the central mass. Most notably, 
Neufeld & Maloney (1995) argue that the outer edge of 
high- velocity emission is due to a transition from a molec- 
ular to an atomic disk at this radius. We discuss this in 
detail in the next section. To date there is no compelling 
explanation for the inner edge of high-velocity emission. 
Here we speculate that the same amplification and absorp- 
tion considerations that make the bowl model attractive 
as an explanation for the systemic geometry play a role in 
the inner cutoff for the high- velocity emission. Figure 21 
shows the disk cross section in the plane of the sky. The 
vertical dashes mark the observed extent of high-velocity 
emission. Water molecules inside of a certain radius will 
have an unabsorbed view of the background jet emission 
and, according to the arguments of Section 4.2, will gener- 
ate strong maser emission perpendicular to our LOS (i.e., 
"their" systemic emission). If these masers are at least 
mildly saturated, as appears to be the case for 'our' sys- 
temic masers (e.g. Modjaz et al. 2005), then competitive 
gain arguments would suggest a corresponding quenching 
of any high-velocity emission along our LOS from within 
this radius. The precise radius at which the background 
jets come into view depends on the orientation and extent 
of the jets. But from Figure 21, which shows the observed 
extent of continuum emission, it seems plausible that there 
is indeed a connection between the purported location of 
'their' systemic emission and 'our' high-velocity emission. 

5.3. X-ray Obscuration and the Accretion Rate 
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For di/dr = 0.034 mas -1 , the accretion disk traced by 
the masers rises in front of and obscures the central engine 
at a disk radius, r x , of 8.3 mas or 0.29 pc (see Figure 10), 
where r x = i' Q j (di/dr) and i' a is the inclination at r = 0. 
The formal fractional uncertainty in this radius is about 
30% and is dominated by the uncertainty in di/dr. How- 
ever, for a disk in which the systemic masers lie along a 
tangent to the disk, r x can be written in terms of r s and 
i s instead of i' Q and di/dr as r x = r s [l + i s /tan(i s )), where 
r s and i s are the radius and inclination (0° for edge-on) 
of the systemic masers. In other words, if we assume that 
the systemic masers are clustered around the bottom of a 
bowl along the near edge of the disk, then the X-ray obscu- 
ration radius is independent of di / dr and depends only on 
other relatively well-constrained parameters of the model. 
We estimate that in the context of the bowl model, r x is 
known to a fractional uncertainty of about 5%. 

From Equation 12, the shielding column density through 
a Shakura-Sunyaev isothermal accretion disk is 

iV 2 4 = 1.3 x 10 -3 cm - 2 . (17) 

Here, N24 is the shielding column in units of 10 24 cm~ 2 , 
M„ 5 = 10 5 M, M 8 = 10~ 8 M o , r pc is the disk radius in pc, 
and c S7 = c s /7kms _1 is the isothermal sound speed. The 
column is computed for an obliquity angle cos -1 ^ with 
respect to the disk normal (fj, = h(r, 9) ■ f) and assuming 
solar abundances. From the bowl model the obliquity an- 
gle is 74 degrees and [i = 0.27. Neufeld & Maloney (1995) 
showed that as a result of the warp in the NGC 4258 disk, 
heating from direct exposure to central X-rays dominates 
viscous heating in the outer disk, and this must be taken 
into consideration in evaluating Equation 17. More pre- 
cisely, they showed that for an externally heated Shakura- 
Sunyaev disk, the midplane pressure scales as r~ 3 . As a 
result there exists a critical radius, r cr , beyond which the 
disk is atomic throughout: 



0.04L 4 i ' 42b (M_ 5 /a) 



0.809.-0.383 



0.617, 



-1.19 



pc, 

where L41 is the 2-10 keV X-ray luminosity in units 
of 10 41 ergs _1 . Numerical simulations by Neufeld, Mal- 
oney, & Conger (1994) indicate that the temperature in 
the predominantly atomic disk beyond r cr rapidly reaches 
an asymptotic value of about 8000 K, corresponding to 
c S7 = 1. Their models also indicate that within r cr , the 
temperature of the dwindling surface atomic layer ranges 
from about 5000 to 8000 K degrees depending on depth 
and proximity to r cr , while the midlplane molecular layer 
temperature ranges from about 400-1000 K. 

Equations 17 and 18 can be used to estimate the X-ray 
shielding column at r x as a function of M_ 5 /a. For each 
accretion rate, Equation 18 is first used to calculate r cr 
along the 6 = 270° radial for c Sr = 1, L 41 = 0.8 (Frus- 
cione et al. 2005), and di/dr — 0.034 mas -1 . If for a given 
Ms/a, r x > r cr , then the disk is entirely atomic at r x 
and N24 is computed from Equation 17 using the asymp- 
totic sound speed of 7kms~ 1 and /j, — fj, rx (fa 0.27). If on 
the other hand r x < r cr , then the X-rays must also transit 
cooler molecular gas in the disk midplane. 

Figure 22 shows N 24 (Ma~ 1 ) at r x = 8.3 mas. The 
abrupt jump in N 24 at M w 10~ 4 a M Q yr _1 marks the ac- 
cretion rate at which r cr passes outside of r x and the disk 



column becomes increasingly dominated by cool, dense 
molecular gas in the disk midplane. We have made no 
effort to model the transition from predominantly molec- 
ular to predominantly atomic gas. However, we note in 
passing that the actual transition is likely to be extremely 
rapid as a result of the strong radial dependence (r -3 ) 
of the disk midplane pressure. Furthermore, Figure 22 is 
based on the assumption that inside of r cr , all the disk 
mass is molecular. While a surface layer of atomic gas 
will in fact persist inside the critical radius, it is reason- 
able to approximate that the column will rapidly become 
doiminated by the dense midplane molecular gas. 

The X-ray column has been observed and estimated nu- 
merous times since 1993 (see Table 2 of Fruscione et al. 
2005 for a summary), and the dotted lines in the figure 
show the range of observed values. Standard AGN models 
suggest that this column is due to a thick dusty torus 1 
to 10 pc from the central engine. With this in mind, we 
treat the observed column as an upper limit on the col- 
umn through the sub-parsec disk. For r cr > r x the disk 
column is 50 — 100 times greater than the largest observed 
column. Hence, for di/dr = 0.034 mas^ 1 , consistency with 
the observed X-ray column requires that the disk be atomic 
at the obscuring radius of 8.3 mas (0.29 pc), translating to 
an upper limit on the accretion rate of <~ 10~ 4 a M yr~ x 
for c S7 = 1 and L41 = 0.8. As discussed above, these 
conclusions are robust to uncertainties in di/dr if the sys- 
temic masers are indeed distributed about the tangent to 
the disk. The upper limit on M assumes a disk filling 
factor of order unity at r x . The unbroken nature of the 
systemic spectrum suggests this is a reasonable approxi- 
mation at <~ 4 mas. Furthermore, the relative stability (to 
factors of a few; Fruscione et al. (2005)) of the observed 
column over the last ten years implies a fairly smooth in- 
tervening medium. Finally, our analysis implicitly treats 
the central X-ray emission as arising from a point source. 
Observed X-ray variability on timcscales of hours (Fiore et 
al. 2001) implies an X-ray source of less than 100 r g corre- 
sponding to an angular size of ~ 10' at the distance of the 
innermost maser. Thus, these conclusions should not be 
sensitive to spatial structure in the central X-ray source. 

Neufeld & Maloney (1995) estimate the accretion rate 
through the NGC 4258 outer disk by postulating that it 
is the transition from a molecular to an atomic disk at r cr 
that defines the radius, r outl of the outermost high- velocity 
mascr. Using Equation 18 and setting r cr = 8.1 mas they 
find M = 7 x 10 _5 aM Q yr _1 . A more conservative ap- 
proach is to require merely that the disk be molecular 
over the entire range of maser emission, leading to a lower 
limit on the accretion rate. Figure 23 shows r cr (Ma _1 ) 
along the midline for the position angle warp discussed in 
the previous section. In order to accommodate both the 
manifest existence of maser emission at r out and the ob- 
served X-ray column at r x , the transition from molecular 
to atomic gas must occur outside of 8.1 mas (0.28 pc) along 
the disk midline and inside of 8.3 mas (0.29 pc) along the 
LOS. For c Sr = 1 and L41 = 0.8, this implies accretion 
rate upper and lower limits of ~ 10~ 4 cn M & yr~ x . Given 
the central mass of 3.79 x 10 7 M©, this translates to an 
accretion rate in Eddington units of ~ 10~ 4 a, assuming 
a canonical efficiency of 0.1. For Lboi ~ 10 42 ergs _1 , we 
estimate a radiative efficiency of > 0.2 for a < 1. 
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The conclusion that the disk transitions from molecu- 
lar to atomic gas between r out and r x is robust in the 
sense that a predominantly molecular disk would lead to 
substantially more column than observed and an atomic 
disk could not support masers. Furthermore, our estimate 
that the transition occurs at ~ 8.2mas= 0.28 pc is rela- 
tively insensitive to details of the fitting process as long 
as it is indeed the case that the systemic masers lie along 
the bottom of a bowl and the high- velocity masers trace a 
position angle warp. The accretion rate estimate, on the 
other hand, is sensitive to estimates for the X-ray lumi- 
nosity and column (both known to be variable by factors 
of a few; Fruscione et al. 2005), the disk sound speed, and, 
more generally, to the applicability of the equations for an 
externally heated thin disk to the outer NGC 4258 disk. 
With this in mind, we refrain from placing formal errors 
on the accretion rate. 

The preceding analysis validates the proposal of Neufcld 
& Maloney (1995) that the high-velocity emission is trun- 
cated at large radii by a transition from molecular to 
atomic gas. In addition, it suggests that a substantial 
fraction, if not all, of the observed X-ray column is due 
to the warped disk at 0.29 pc (Figure 22). The standard 
model of low-luminosity AGN (Antonucci & Miller 1985) 
asserts that Seyfert Is and 2s are the same class of object 
observed from different perspectives. More recent work 
suggests that HBLR and non-HBLR Seyfert 2s may in fact 
be distinct classes of objects, with the HBLR Seyfert 2s 
being the misaligned counterparts to the face-on Seyfert Is 
(c.f. Veilleux 2003). While it is generally accepted that at 
least some Seyfert 2s are Seyfert Is shrouded in intervening 
gas, the precise nature of the absorbing material is still de- 
bated. Broad emission-line region clouds, sub-parsec-scale 
warped disks, 1-10 pc thick molecular tori, and larger-scale 
wind structures have all been considered. Recent studies 
of column variability in Seyfert 2 galaxies suggest inter- 
vening material quite close to the central engine (Risaliti, 
Elvis, & Nicastro 2002). NGC 4258 is just one object, and 
as such generalizations cannot be drawn from it. Never- 
theless, our analysis, together with the X-ray observations 
of Fruscione et al. (2005), presents compelling evidence 
for at least one HBLR Seyfert 2 in which most, if not all, 
of the X-ray obscuration arises in a subpar sec- scale warped 
disk. 

If the X-ray column indeed arises in the disk, then 
the expected timescale, t var , for variability in the X- 
ray column is H/v rot , where we have associated the disk 
thickness with the characteristic scale of the intervening 
material. For an isothermal disk, H = rc s /v rot , and 
t var = rcs/vrot 2 . For the NGC 4258 disk at 0.29 pc, 
v rot ~ 750 km s -1 and, taking c s = 7kms _1 , we estimate 
t var fa 3.6 yrs. The observed column variability of ~ 60% 
over ~ 5 months (Fruscione et al. 2005) is approximately 
consistent with this prediction, but may indicate inhomo- 
geneities on a scale somewhat less than the disk thickness. 

The existence of water masers requires 300 < T m < 
1000 K and 10 s < p m < 10 10 cm" 3 in the NGC 4258 disk, 
where T m and p m are the temperature and density of the 
masing gas. This temperature is consistent with the mod- 
els of Neufeld et al. (1994) for an externally irradiated thin 
disk. The density of the maser layer is estimated by as- 
suming that the masers lie approximately at the midplanc 



of the disk. From Equation 12 the midplane number den- 
sity, n H2 = p m id/m H2 can be written 

n H2 = 3700^^- cm" 3 . (19) 
r i c s i a 

For M = 1 x 10- 4 a M Q yr" 1 and c s = 2.3 kms" 1 ^ = 
800 K), nH2 ranges from 7 x 10 8 cm -3 at the inner radius 
of 0.14 pc to 0.9 x 10 8 cnT 3 at 0.28 pc, the radius of the 
outermost high-velocity feature. Hence, the accretion rate 
derived from the disk geometry and X-ray column provides 
conditions suitable for water masers in the disk midplanc 
over the entire envelope of observed maser emission. 

The same disk material that obscures the central X-ray 
emission will also obscure any radio emission arising in the 
central engine. Herrnstein et al. (1998) report a 220 ply 
upper limit on 22 GHz emission coincident with the cen- 
tral mass, and this limit has been used to place constraints 
on the accretion mode and rate through the disk (Herrn- 
stein et al. 1998; Gammie, Narayan, & Blandford 1999; 
Yuan et al. 2002). The 22-GHz optical depth due to free- 
free absorption in the disk at r x will be 1-2 (Herrnstein, 
Greenhill, & Moran 1996), and hence we speculate that 
the actual upper limit on central-engine 22 GHz is likely 
1-2 times larger than the observed 220 pJy. 

Finally, we note that any association between the ac- 
cretion rate in the maser disk and that in the NGC 4258 
central engine implicitly assumes that mass loss via winds 
can be ingored and that the disk is in a steady state, 
such that dM{r)/dr = 0. Gammie et al. (1999) argue 
that in the absence of external torques, steady state is 
achieved after a viscous timescale given approximately by 
t v ~ 0.26a- 1 (200 K/T)(r/0.2 pc) 1 / 2 Gyr. For T = T m = 
800 K and a > 0.3, t v < 0.2 Gyr at the radii of the masers. 

5.4. The NGC 4258 Central Engine 

Ho (2003) lists the defining properties of Low luminosity 
AGN. LLAGN are highly sub-Eddington, and their spec- 
tral energy distributions (SEDs) lack the 'big blue bump' 
of more powerful AGN. Instead they typically show an ex- 
cess of IR emission, perhaps even peaking in this region 
of the spectrum. They tend to be relatively strong radio 
emitters, and they possess narrow rather than relativis- 
tically broadened iron ka emission. NGC 4258 satisfies 
all of these criteria, especially if we accept that the steep 
(s = 1.1±0.2) optical/UV spectrum of Wilkes et al. (1995) 
is in fact an extension of the steep (s = 1.4 ±0.1) IR spec- 
trum of Chary et al. (2000), and that the SED therefore 
peaks in the IR. 

The generic LLAGN SED as described above shows in- 
dications that at least three distinct emission mechanisms 
may be at work. First, the IR excess (as well as nu- 
merous examples of double-peaked emission lines; Ho et 
al. 2000) is consistent with black body emission from a 
standard optically thick, geometrically thin, radiatively 
efficient Shakura-Sunyaev disk. However, the absence of 
an optical/UV peak and of relativistically broadened ka 
lines indicates this disk does not extend too close to the 
black hole. Next, the representative LLAGN SED seems 
broadly consistent with emission from a class of hot, quasi- 
spherical, inefficient accretion flows, known collectively as 
radiatively inefficient accretion flows (RIAFs; see Quaetert 
2003 for a recent summary). Accretion of this type natu- 
rally explains the low bolometric luminosities and inferred 
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low accretion efficiencies of LLAGN. Finally, the character- 
istic radio excess is most readily explained by compact jet 
emission at the central engine (Falcke & Biermann 1995, 
1999; Anderson et al. 2004). 

Hence, a plausible model for the LLAGN central en- 
gine (c.f. Ho 2003; Qataert 2003) has an outer thin accre- 
tion disk truncated at <~ 100 — 1000r ff by a hot, inefficient 
RIAF, which in turn channels some fraction of the accreted 
energy into a compact jet. While it is perhaps uncontro- 
vertial that one or all of these mechanisms play a role 
in most LLAGN central engines, the neat itemization of 
the previous paragraph is probably misleading. The RIAF 
models can be enhanced to accommodate the radio excess, 
the compact jet models can be extended to accommodate 
essentially the entire SED, and the thin disk can be left 
out altogether. 

Several attempts have been made to model the 
NGC 4258 spectrum. Gammic et al. (1999) propose an 
accretion rate of ~ O.OlMQyr -1 through an advection 
dominated accretion flow (ADAF) that transitions to a 
thin disk at 10 — 100 r g . The model relies on relatively 
high accretion rates through an outer thin disk to explain 
the observed IR luminosity and the extreme inefficiency 
of the central ADAF to accommodate the significant sub- 
Eddington luminosity of the system. However, this model 
appears to be inconsistent with the subsequent detection 
of a steep (s = 1.4 ± 0.1) IR spectrum by Chary et al. 
(2000). Yuan et al. (2002), on the other hand, propose 
that essentially the entire NGC 4258 LLAGN spectrum 
arises in a compact jet. They demonstrate broad agree- 
ment with the observed SED, and in particluar the steep 
IR spectrum, for jet mass loss rates between 2.9 x 10~ 4 
and 7 x 10~ 5 M Q yr _1 . They speculate that disk accretion 
rates of 10~ 3 — 10~ 4 M Q yr _1 are needed to feed the jets, 
and show that for such low accretion rates, an outer thin 
disk and central ADAF contribute negligably to the jet 
sprectrum. In this jet-dominated model, the radiatively 
inefficient ADAF serves primarily as a means to channel 
large amounts of accretion energy to the radiatively effi- 
cient jets. Our estimate for the NGC 4258 accretion rate 
is in reasonable agreement with the jet-dominated models 
of Yuan et al. (2002). In NGC 4258 we know that all the 
detected 22 GHz emission arises in compact jets and ap- 
parently not in a central RIAF. Furthermore, we appear 
to be able to place relatively tight constraints on the ac- 
cretion rate at 0.28 pc, and this may aid in untangling the 
relative contributions of jet and RIAF. 

6. CONCLUSIONS 

In this paper, we have proposed a geometry for the sub- 
parsec disk in NGC 4258 based primarily on its ability to 
explain a large number of oherwise puzzling features of the 
system. In particular, an inclination warp explains: 

• The sub-Keplerian projected rotation curve of the 



high- velocity masers and implies a purely Keplerian 
underlying disk. 

• The narrow declination spread of the systemic 
masers, and their apparent tight clustering around 
a radius of ~ 4 mas, corresponding to the bottom 
of a bowl along the near edge of the disk. 

• An apparently recurring flare in the systemic 
masers <~ 20kms _1 red of the systemic velocity of 
the galaxy. 

• In conjunction with observed X-ray column, the 
outer cutoff of high-velocity emission. 

• By analogy with amplification and obscuration ar- 
guments for the systemic masers, the inner cutoff 
of the high-velocity emission. 

While the model is compelling in its ability to explain a 
variety of aspects of the system, we have not demonstrated 
that it is the best model in the formal maximum likelihood 
sense. In particular, we have shown that there remains a 
2cr tension between the disk center predicted purely from 
the high-velocity masers, and the center derived from the 
high- velocity masers and the flare location. The resolution 
of this tension requires global fits, perhaps incorporating 
additional details not considered here (such as non-circular 
orbits). 

Finally, the inclination-warp model makes several pre- 
dictions, most obviously that systemic maser emission can- 
not arise significantly further south of disk center than 
0.55 mas, and that all high- velocity emission must occur 
within ~ 8.1 mas. This has held true up to this point, and 
will be easily verifiable in upcoming datasets. In addition, 
it predicts that future flares should occur ~ 20kms~ 1 to 
the red of the systemic velocity. If it is indeed the case that 
the inner edge of high- velocity emission is due to compet- 
itive gain with stronger, amplified emission in the plane 
of the sky, then this is not a hard edge like the outer cut- 
off. To the extent that the systemic spectrum is variable 
and occassionally weak at certain velocities, we might ex- 
pect to see faint high- velocity emission inside of the ususal 
<~ 0.16 pc inner cutoff. Finally, we point out that in this 
model there is no a priori reason that the systemic masers 
must lie along the bowl bottom. Rather, this is where 
such emission will be brightest. We speculate that there 
could be additional much weaker low- velocity emission dis- 
placed from the bowl bottom. The most likely place to find 
such emission would be <~ 20kms~ 1 to the red of the sys- 
temic velocity, where the velocity gradient along the LOS 
is smallest. 

We thank R. Herrnstein for constructive discussions and 
for her tireless help in bringing this project to completion. 
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Table 1 

NGC 4258 VLBI Observing Log 



Experiment 


Observation 


Antennas*- 1 ^ 


Frequency 


Synthesized 


Sensitivity^ 


Code 


Date 




Setups) 


Beam' 3 ) 




BM19 


19940426 


V+Y 


L 


1.1x0.6 @ 9° 


6.0 


BM36a 


19950108 


V+Y 


L 


0.8x0.7 @ -15° 


6.5 


BM36b 


19950529 


V+Y 


L 


1.0x0.7 @ -10° 


5.2 


BM56b 


19960921 


V+Y+G+E 


L 


1.3x1.1 @ 20° 


3.0 


BH25a 


19970307 


V+Y+G+E 


C 


0.7x0.3 @ 15° 


3.0 


BH25b 


19970323 


V+Y+G+E 


C 


0.6x0.3 @ 15° 


2.3 


BH25c 


19970407 


V+Y+G+E 


c 


0.6x0.3 @ 15° 


2.6 



1 V: VLBA; Y: phased VLA; G: NRAO 140-foot; E: Effelsberg 100-meter. 
2 : Line mode; C: Continuum mode (see text). 
3 mas x mas © position angle. 
4 mjy/km/s/beam. 
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Table 2 
Rotation Curve Models 



Model 


Relativity? 


X a 


vo 


M 


Xr 


Parameter Value 


Av b 






(mas) 


(km s _1 ) 


(xlO 7 Mq) 






(km s- 1 ) 


Keplcrian 


No 


-0.29 ±0.03 


488 ±2 


3.85 ±0.01 


1.24 






Keplerian 


Yes 


-0.33 ±0.02 


486 ±2 


3.84 ±0.01 


1.24 






Power Law 


Yes 


-0.19 ±0.04 


475 ±3 


3.56 ±0.01 


1.09 


0.48 ±0.01 c 


9.0 


Central Cluster 


Yes 


-0.19 ±0.05 


475 ±4 


3.74 ±0.02 


1.07 


5.7(±1.4) x 10 11 d 


10.9 


Massive Disk 


Yes 


-0.20 ±0.04 


476 ±3 


3.80 ±0.01 


1.10 


6.5(±1.3) x 10~ 2 c 


10.1 


Inclination Warp 


Yes 


-0.24 ±0.03 


479 ±3 


3.78 ±0.01 


1.08 


0.04±0.01 f 


9.6 


Inclination Warp 8 


Yes 


-0.191 


473.5 


3.794 


1.19 


0.034 f 


9.64 



"Offset with respect to reference maser at 510 km s -1 LSR. 

'Maximum deviation from Keplerian rotation between innermost and outermost masers, positive representing a flattening in the projected 
rotation curve. 

C £ (Vrot oc ) 

V (M s pc- 3 ) 

e M/a (M Q yr- 1 ) 

fdi/dr (mas -1 ) 

9 di/dr bootstrapped global model 
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LOS Velocity (km s-1 ) 

Fig. 1. — Top: Warped-disk model with masers and continuum superposed. Bottom: Total power spectrum of the NGC 4258 maser, with 
best-fitting Keplerian rotation curve. Maser data are from VLBA epoch BM36b, continuum data are from epochs BH25a-c (see Table 1). 
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o.o 



500 



3.75 3.80 3.85 

M„ (1 7 Solar Masses) 

Fig. 2. — Probability density functions for each of the disk center parameters. The solid lines assume a pure Kcplcrian rotation in a constant 
inclination disk, while the dashed lines are derived for a model incorporating an inclination warp. The shaded region in the middle plot shows 
the 1(7 confidence interval for the galaxy systemic velocity, as measured from optical observations. 



1050 p 




4567845678 
Radius (mas) Radius (mas) 

Fig. 3. — Best-fitting high- velocity rotation curve models derived using epochs 1—4 in Tabic 1 (solid lines). Open circles mark positions 
and velocities of red maser features, filled squares of the blue features. Residuals are shown in the lower panels. Panel a shows a Keplerian 
model, neglecting relativistic effects. All subsequent panels incorporate relativity, b) Keplerian model, c) Power law model (v ro t oc r - '). d) 
Central cluster model, e) Massive disk model, f) Inclination warp model. 
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-4 -3 -2 -1 -0.05 -0.02 0.01 0.04 0.07 0.10 

log[M/a (M s yr~')] di/dr (mas - ') 



Fig. 4. — Relative x 2 f° r f° ur different refinements to a Keplcrian, constant-inclination disk. Clockwise from upper left they are: 1) 
Power law rotation parameterized as v oc r~^; 2) a central cluster with core density p c ; 3) a massive disk, parameterized in terms of the 
accretion rate M; and 4) a disk with a radial gradient in inclination of di/dr. The solid curves show the effect of these perturbations on a 
constant-inclination disk, the dashed and dotted curves on a disk incorporating di/dr = 0.04 mas -1 and di/dr = 0.034 mas -1 , respectively. 




-3.0 -2.5 -2.0 -1.5 -1.0 -0.5 0.0 -0.02 0.00 0.02 0.04 0.06 0.08 

log[M/a (M Q yr -1 )] di/dr (mas -1 ) 

Fig. 5. — Probability density functions for each of the rotation curve perturbations considered in Figure 4. For the power law model 
(upper-left panel), the solid curve shows the distribution of £ for a constant inclination disk, while the dashed and dotted curves show the 
distribution for a disk incorporating di/dr = 0.04 mas -1 and di/dr = 0.034 mas -1 , respectively. 



o r 




Fig. 6. — Simultaneous relative \ 2 intervals for M and di/dr. Contour steps are 1, 1.5, 2, 3, 4, 5, 7, 11, and 15. The strong correlation 
between the effects of a massive disk and an inclination warp on the rotation curve is evident. 
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-0.1 -0.2 
X Offset (mas) 

Fig. 7. — Relative RA versus declination for the systemic masers for epochs BH25a (+), BH25b (*), and BH25c (diamonds). The best-fitting 
line corresponds to a position angle of 5.6°. The error bars include a non-instrumental noise of 8 /ias added in quadrature to achieve a reduced 
X 2 
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Z Offset (mas) 



Fig. 8. — Cross-section of a constant-inclination disk. The LOS is along the —z direction. The northern shaded oval marks the observed 
location of northern jet emission, and the southern oval is the reflection of this through the disk center. The asterisk marks the observed 
location of the southern jet emission. The small vertical dash at a Y Offset of —0.55 denotes the position of the systemic masers in the disk, 
drawn to scale in Y. It is difficult to reconcile the narrow Y confinement of the systemic masers with background amplification models and 
the large radial spread of the high- velocity masers. 
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Fig. 9. — Relative RA versus LOS velocity for the systemic masers for epochs BH25a (+), BH25b (*), and BH25c (diamonds). The fits 
were performed assuming x = —0.19 mas and </> = 5.6° (see text for details), a) and b) Linear fit assuming cr r /r = 0.015, yielding a slope of 
270 km s -1 mas -1 and a reduced \ 2 of 40. c) and d) Fit using oy/r = 0.13 giving a slope of 266 km s -1 mas -1 and a reduced \ of 1. The 
insets show in greater detail the portion of the disk that dominates the \ 2 - 




> 



10 5 -5 -10 

Z Offset (mas) 

Fig. 10. — Cross-section of an inclination-warped disk with di/dr = 0.034 mas -1 . The inclination angle measured downward from the LOS 
is i = 0.28 + 0.342, where z is in milliarcscconds and i is in radians, and the shape of the warp is y = z sini. The LOS is along the — z direction. 
The northern shaded oval marks the observed location of northern jet emission, and the southern oval is the reflection of this through the 
disk center. The asterisk marks the observed location of the southern jet emission. The small vertical dash at a Y offset of —0.55 mas denotes 
the position of the systemic masers in the disk, drawn to scale in Y. 




Fig. 11. — Schematic drawing of the bowl bottom. Both the thickness of the bowl as well as the curvature have been exaggerated. The 
thick line on the upper surface of the disk represents a thin layer of gas ionized by x-ray emission from the central engine. The shaded region 
demarcates the expected region of brightest maser emission. See text for a discussion of points A, B, and C. 
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Fig. 12. — Top: Relative RA versus declination for the systemic masers for epochs BH25a (+), BH25b (*), and BH25c (diamonds). Error 
bars are the instrumental uncertainties. The solid line marks the bottom of the bowl for di/dr = 0.034 mas -1 . The dotted line corresponds 
to a la maser layer thickness of 15 fias. Bottom: Radius as a function of x for the inner edge A, bottom B, and outer edge C of the bowl 
(as identified in Figure 11) assuming a maser layer thickness of 15/ias. 




0.20 0.10 0.00 -0.10 -0.20 -0.30 -0.40 -0.50 

X Offset (mas) 

Fig. 13. — LOS velocity versus x offset for the systemic masers for epochs BH25a (+), BH25b (*), and BH25c (diamonds). The solid lines 
show model LOS velocity as a function of x for the positions indicated in Figure 11, assuming a maser-layer thickness of 15 /ias, and using 
v = 473.5kms~ 1 and x = — 0.19mas. The arrow indicates the maser features associated with the flare, detected in each of the three epochs. 
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X Offset (mas) 



Fig. 14. — Contours of constant LOS velocity relative to v = 473.5 km s —1 along the bottom of the bowl. The LOS is along the —z direction, 
and hence vertical contours indicate a zero velocity gradient across the bowl along the LOS. The model used x„ = —0.29 mas. The roughly 
horizontal line marks the position of the bottom of the bowl (B in Figure 11). The intersection of this line with the zero contour locates the 
projected near edge of the disk (x' a ) at ~ —0.15 mas. 
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Fig. 15. — Time series of Haystack single-dish spectra of systemic maser emission. Labels indicate the epoch of observation (yyddd). A 
vertical line marks approximately the velocity of the flaring spectral feature at the end of the monitoring, ~ 2 kms -1 above its initial velocity. 
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Fig. 17. — Systemic spectra from VLBA epochs BH25a, BH25b, and BH25c showing the 1997 flare. These were observed on days 66, 82, 
and 97 of 1997, respectively. 
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Fig. 18. — Systemic spectra of the 1996 and 1997 flares. The solid line shows the 1996 flare as detected by the VLA on day 89 of 1996 
(from Bragg et al. 2000). The other two spectra show single-dish detections of the 1997 flare at the start (dotted) and end (dashed) of its 
light curve, on days 77 and 175, respectively, of 1997. 
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Fig. 19. — Model disk center parameters for a variety of fits. The diamonds mark the disk center from fits incorporating only the high- 
velocity data: open for a purely Keplerian model, filled for an inclination-warp model. The squares are derived from fits incorporating both 
the high-velocity and systemic data. The open square is from the global fit of Herrnstcin (1997). The solid square is derived from the 
flare-constrained model of Section 4.3. The solid line marks the locus of disk centers for which the velocity convergence point lies along the 
observed systemic emission. The dotted lines demarcate the la confidence interval for the galaxy systemic velocity as measured from optical 
observations. 
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Fig. 20. — Three perspectives of the best-fitting disk incorporating both a position angle warp and an inclination warp with di/dr = 
0.034 mas -1 . The line of sight is shown extending beyond the outer edge of the disk. The disk midline, the = 0° and 9 = 270° radials, 
and the locus of disk tangents are overlaid on the disk. The masers, from all the epochs in Table 1, have been positioned according to their 
observed x and y, and assuming that the high-velocity features lie along the midline and the systemic features lie along the bowl bottom. 
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Fig. 21. — Cross-section of the disk in the plane of the sky, along the disk midline, showing the position angle warp. The warp has been fit 
using VLBA epochs BM19, BM36a, BM36b, and BM56b, and is given by <j> = 0.425 - 0.088r + 0.0038r 2 . The small vertical dashes mark the 
envelope of observed high-velocity maser emission. The northern shaded oval marks the observed location of northern jet emission, and the 
southern oval is the reflection of this through the disk center. The asterisk marks the observed location of the southern jet emission. 
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Fig. 22. — X-ray column density through the disk as a function of MaT x for L41 = 0.8. The column is calculated at a radius, r x , of 8.3 mas 
(0.29 pc), at which point the disk rises in front of the central X-ray source for di/dr = 0.034 mas -1 . For M < lO~ 4 aM0yr _1 , the disk at r x 
is atomic throughout. As the accretion rate passes through lO _4 oM0yr _1 the disk at r x transitions to cool, molecular gas, and this leads 
to the sudden jump in the column. The two dotted lines mark the range of observed X-ray columns in the XMM data of Fruscione et al. 
(2005), and provide an upper limit to the disk column. 
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Fig. 23. — Atomic-to-molccular phase transition radius along the disk midline as a function of MoT 1 for c s? = 1 and L41 = 0.8. The 
dashed line corresponds to the radius of the outermost high- velocity maser, and provides a lower limit to the actual phase transition radius. 



